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1.  INTRODUCTION 


This  report  summarizes  the  five  years  of  AFRL  supported  research  on  the  characteristics 
and  structure  of  the  ionosphere  as  it  impacts  the  operational  requirements  of  the  US  Air  Force. 
Three  major  areas  of  our  research:  equatorial  scintillation,  OpSEND  and  CalibrationA^ alidation 
of  DMSP  UV  measurements  are  addressed  in  this  report.  We  also  include  a  description  of  the 
extensive  support  effort  directed  towards  the  operational  facilities  of  the  Air  Force. 

2.  EQUATORIAL  SCINTILLATION 

2.1  Introduction 

It  is  well  established  that  on  many  days,  after  sunset,  the  equatorial  F-region  becomes  unstable, 
resulting  in  amplitude  and  phase  disturbance  of  trans-ionospheric  satellite  signals.  The  principal 
problem  here  is  to  determine  which  days  will  be  effected  and  which  will  remain  undisturbed. 
This  prediction  problem  involves  the  understanding  of  the  trigger  mechanism  and  the  conditions 
necessary  for  the  rapid  growth  of  these  Rayleigh-Taylor  (R-T)  instabilities.  This  project 
involves  analysis  of  scintillation  phenomena  using  one  year’s  satellite  scintillation  data  from 
Ancon,  Peru  (Jicamarca).  Two  geostationary  satellites  were  used.  These  satellites  transmit  at  a 
fi-equency  of  about  250  MHz.  One  is  over  the  Pacific  Ocean  to  the  west  of  Ancon  and  another  is 
over  the  Atlantic  Ocean  to  the  east.  The  intersection  of  the  line-of-sight  path  fi’om  the  satellites 
to  Ancon  through  the  F-region  is  shown  in  Figure  1.  Figure  1  also  shows  for  a  second 


Longitude 

Figure  1 .  Location  of  the  Ionospheric  Sub-Satellite  Points  from  the  Atlantic  and  Pacific 
Satellites  are  shown  as  a  function  of  altitude  fi-om  300  to  900  km.  The  near  vertical  curves  show 
the  magnetic  flux  tubes  aligned  along  the  local  magnetic  meridian  in  this  region. 
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receiving  station  at  Antofagasto,  Peru,  near  the  southern  anomaly  region,  the  ionosphOTC 
intersection  points  from  the  same  two  satellites.  The  indicated  points  show  the  intersection 
location  as  a  function  of  height,  from  300  to  900  km  altitude. 

The  system  at  Ancon  continuously  records  each  satellite  signal,  and  the  scintillation 
index  S4  was  used  to  indicate  the  presence  of  scintillations.  Scintillations  typically  begin  within 
one  hour  after  sxmset  and  can  continue  for  many  hours  during  the  night.  Along  with  the 
scintillation  measmements,  the  DPS  ionospheric  soimder  at  Jicamarca  was  operated  on  a  one  half 
hour  schedule  during  the  year.  These  routine  sounder  measurements  and  ionogram  inversions 
made  it  possible  to  follow  density  and  height  variations  that  are  associated  with  the  onset  of 
scintillation.  After  sunset,  spread-F  on  the  ionograms  appears  as  the  dominant  feature  associated 
with  scintillation.  This  investigation  attempts  to  develop  the  relationship  between  the  early 
stages  of  spread-F  and  the  development  of  full-blown  scintillation. 

2.2  Scintillation  Statistics 

Scintillation  begins  after  F-layer  (300  km)  sunset,  typically  from  30  minutes  to  90 
minutes  after  the  F-layer  is  no  longer  illuminated.  Figure  2  combines  all  the  data  from  the  two 
satellites  to  show  the  onset  time  of  scintillation  (S4  >  0.6)  relative  to  the  F-layer  sunset  for  the 
entire  year  of  1998. 
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Figure  2.  Occurrence  of  scintillation  relative  to  sunset  at  Ancon  using  geostationary  satellites  for 

1998. 
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Two  factors  become  obvious  from  these  data.  First,  between  day  100  and  240 
^proximately  Apnl  through  August)  the  occurrence  of  strong  scintillation  is  relatively  rare 
This  IS  a  well-known  characteristic,  though  the  period  of  minimum  scintillation  occurrence 
anes  \^th  geographic  location,  apparently  depending  on  the  local  magnetic  declination  The 
second  factor  is  best  examined  using  the  histogram  (Figure  3)  of  the  onset  times  for  all  these 
data,  excluding  days  1 00  through  240. 


Figure  3.  Histogram  of  the  onset  time  relative  to  sunset  for  strong  scintillation  at  Ancon  for  the 

year  1998  (all  days  except  100  through  240). 

mA/i-  ni^'ority  of  Aese  data  lie  between  0  and  60  minutes  after  F-layer  sunset  with  a 
median  value,  bas^  on  the  entire  distribution,  of  48  minutes.  Time  delays  a  long  as  three  hours 
were  obs^ed  and  less  than  five  percent  of  the  events  occurred  before  sunset.  As  will  be 
discussed  later  in  this  report,  the  growth  factor  for  the  R-T  instabilities  depends  primarily  on  the 
chai^ter  of  fee  F-region  lomzation  near  fee  bottom  of  fee  layer.  The  first  look  at  fee  chLges  in 
the  height  of  fee  bottom  of  fee  layer  was  shown  to  give  a  good  prediction  of  fee  onset  of 
scmtillation  vnth  a  lead-time  of  one  to  two  hours.  This  was  based^on  fee  enhanced  verticaf 
movement  of  fee  F-layer  feat  begins  just  before  sunset.  This  movement,  an  ExB  vertical  drift 
IS  ^im^ted  to  changes  m  fee  zonal  electric  field.  Figure  4  shows  fee  virtual  height 
h^ges  of  fee  F-layer  for  fee  months  of  October,  November,  and  December  1998.  These  resifes 
followed  the  hei^t  vanations  usmg  a  sounding  frequency  of  5  MHz.  It  was  important  to  use  a 
frequency  tlwt  reflects  from  fee  F-region  during  fee  entire  day.  Although  there  appears  to  be 
the^SS  changes  approximately  one  hour  before  sunset  when 

this  1  T'  varies  considerably  during 

of  to  of  hei^t  increase,  fee  F-layer  maximum  height  reaching  a  virtual  height 

of  500  to  600  km  on  some  mghts  and  changing  only  a  small  amount  on  others. 
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Figure  4.  Jicamarea  virtual  height  variation  at  5  MHz  for  Oct.,  Nov.  and  Dec.,  1998 

Stahshcally,  data  c^  —  " 

ZmLl  alh^de  for  a  "  lI*-Zchoic‘e  SS 

the  orobability  of  scintillation  is  0.20,  while  for  Zpeak  >  400  km  me  Pi^oo^m  y 
"““te  sm.is.ic  is  dm.  if  spread-F  is  ob^ed 

:Sa.t.“s  “of  S’“^cdve  vairsTnce  d,e  tinie  difference  between  d.e« 
occurrences  is  typically  less  than  30  minutes. 


4 


0.8 


...  0.2  •- 

z<400km  •  ■■ 

S4  <  0.4 

0.1  - 


0.0  ■ 


'  r 

300  320  340  360  380  400  420  440  460  480  500  520 

Peak  altitude  after  sunset  (km) 


Figure  5.  Distribution  of  nightly  S4  events  as  a  function  of  peak  F-layer  virtual  height. 

2.3  Instability  Motions 

These  R-T  instabilities  drift  eastward  after  they  are  formed  with  the  ExB  drift  associated 
with  the  downward  electric  field  at  an  average  speed  of  100  m/s.  Analysis  of  many  Jicamarca 
ionograms  indicated  that  the  sounder  could  detect  the  large  bubbles  before  they  reached  the 
overhead  position  and,  in  a  sense,  the  bubbles  could  be  tracked  as  they  approached  the  sounder 
location.  Unfortunately,  the  Jicamarca  sotmder  ran  on  a  30-minute  schedule  that  limits  its  ability 
to  precisely  determine  the  drift  speed.  A  sequence  of  ionograms  is  shown  in  Figure  6,  beginning 
at  00  UT  (F-layer  sunset).  The  bubble  instability  reaches  the  overhead  position  at  0200  UT. 
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Figure  6.  Jicamarca  ionograms  sequence  from  24  February  1997  showing  the  approach  of 

the  bubble  region. 

Combining  the  sounder  data,  the  Julia  radar  and  the  scintillation  index  give  a  composite 
picture  of  the  equatorial  phenomenology  during  this  sunset  period,  as  shown  in  Figure  7.  The 
Julia  radar  at  the  Jicamarca  radar  facility  is  a  subset  of  the  incoherent  scatter  radar  that  operates 
at  lower  power.  At  this  level  it  is  able  to  detect  only  the  stronger  coherent  scatter  components  of 
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Figure  7.  Composite  of  ionosonde  true  height  contours,  the  Julia  radar  and  scintillation 
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irregularity  structure  passing  over  Jicamarca.  However,  the  lower  power  requirement  allows  the 
Julia  radar  to  operate  on  a  more  frequent  basis  and  is  used  here  to  compare  with  the  digisonde 
data.  Superimposed  on  the  plasma  frequency  contours  derived  from  the  ionogram  inversions  for 
the  24-hour  period  are  the  Julia  radar  data  and  plotted  above  is  the  Ancon  scintillation  index  for 
satellite.  On  the  basis  of  these  data,  the  radar  shows  that  bubbles  originates  near  the  bottom  of 
the  layer,  specifically,  very  close  to  the  to  the  region  of  maximum  electron  density  gradient  as 
shown  by  the  diamond  symbol  trace.  This  is  associated  with  the  R-T  instability  growth  factor  as 
discussed  in  the  next  section.  The  Julia  radar  scatter  grows  to  heights  above  the  peak  of  the  F- 
layer.  In  this  example,  the  bubble  structures  arrive  over  Jicamarca  at  around  0200  UT,  with 
bottom-side  spread-F  beginning  approximately  one  hour  earlier.  The  scintillation  appears  to 
begin  with  the  appearance  of  bottom-side  spread-F  but  it  must  be  remembered  that  the 
ionospheric  intersection  point  of  the  scintillation  lies  to  the  west  of  Jicamarca  and  more  closely 
coincides  with  the  bubble’s  passage  before  it  reaches  the  overhead  position. 

Figure  8  illustrates  the  relationship  between  the  bubble’s  eastward  drift  and  the  several 
measuring  techniques  used  here.  Three-dimensional  numerical  ray  tracing  was  used  to 
rmderstand  how  the  sounder  is  able  to  image  the  approaching  bubbles  before  they  reach  the 
overhead  position.  The  model  of  the  bubble  (a  colunmar  region  depleted  of  90  percent  of  the 
ambient  ionization  and  reaching  from  the  bottom  of  the  F-layer  through  the  peak  of  the  layer) 
was  superimposed  on  the  PRISM  ionospheric  model.  The  ray  tracing  was  carried  out  at  7  MHz 
for  several  positions  of  the  bubble  region  with  respect  to  the  soimder  location  (Jicamarca). 

These  calculations  show  the  existence  of  many  points  inside  each  bubble  where  the  ray  is 
orthogonal  with  the  local  magnetic  field  indicating  the  potential  for  scatter  back  to  the  soimder. 
On  the  basis  of  these  calculations  it  expected  that  the  sounder  is  just  able  to  detect  the 
approaching  bubbles  at  a  maximum  range  of  500  km  and  more  easily  as  the  bubbles  get  closer  to 
Jicamarca. 


Figure  8.  Illustration  of  the  relationship  between  eastward  drifting  but  ;s  and  the  several 

measurement  techniques. 
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Finally,  for  several  nights  the  tracking  of  the  bubbles  by  the  sounder  was  used  to 
determine  the  eastward  drift  velocity.  These  measurements  are  shown  in  Figure  9,  plotting 
location  against  time. 


Observed  bubble  speed  using  the  Digisonde 


Figure  9.  Observed  distance  of  equatorial  bubbles  from  Jicamarca  as  a  function  of  time  on 

four  nights. 


The  average  speed  varied  from  80  to  140  m/s  on  these  days. 

Following  Kelley  [1989]  the  R-T  instability  ^owth  factor  is  expressed  below  in  terms  of 
the  zonal  electric  field  intensity  Ex,  the  geomagnetic  field  B,  the  acceleration  of  gravity  g,  the 
ion-neutral  collision  frequency  Vm  and  the  layer  tilt  angle  x. 


/rt 

/RT  ^ 


— cos(r)  +  cos(r)  +  sin(r) 

B  V.  B 


The  length  scale  L  is  given  in  terms  of  the  electron  density  profile  as  i  = 


1  dn^ 
n  dz 


(1) 

(2) 


The  largest  instability  growth  rate  occurs  when  L  is  minimum  or  where  dn/dz  is  largest. 
In  general,  this  occurs  near  the  bottom  of  the  F-layer.  Future  research  involves  determining  from 
each  measured  electron  density  profile  where  the  minimum  L  occurs  and  the  magnitude  of  i  at 
that  point.  For  simplicity,  if  the  layer  tilt  angle  is  assumed  to  be  zero,  then  the  third  term, 
containing  sin  (x),  is  zero  and  the  resultant  expression  becomes: 


E  g 

^4.-^ 

B 


in  J 


(3) 
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The  remaining  important  factors  are  the  magnitude  of  the  zonal  electric  field  and  the  ion- 
neutral  collision  fi-equency  that  decreases  with  increasing  height  of  the  F-layer.  This  explains  the 
increasing  probability  of  scintillation  occurrence  as  the  layer  is  driven  upward  around  sunset. 
The  zonal  electric  field  then  plays  an  important  dual  role.  First,  directly  in  the  equation  of  the 
growth  rate  and  second,  in  driving  the  layer  upward,  thereby  decreasing  the  ion-neutral  collision 
fi-equency.  Both  effects  work  to  increase  the  instability  growth  rate.  The  measurement  of  the 
zonal  electric  field  (vertical  ExB  drift)  is  discussed  in  Section  4  and  the  ion-neutral  collision 
fi’equency  is  available  fi'om  the  sounder  electron  density  measurements  and  the  MSIS  neutral 
atmospheric  model. 

2.4  Drift  Analysis 


On  the  basis  of  the  above  analysis  on  the  R-T  instability  growth  rate,  the  most  important 
factor  is  the  upward  drift  of  the  layer  beginning  just  before  sunset.  The  Digisonde  in  the  “drift 
mode”  is  able  to  measure  a  vertical  drift  velocity.  Radio  sounding  measurements  of  vertical  drift 
has  two  components.  The  first  is  the  motion  of  the  reflection  point,  which  depends  on  the 
sounding  jfrequency,  as  the  layer  moves  upward  (or  downward)  and  the  second  component  is  the 
time  dependent  change  in  the  electron  content  below  the  reflection  level  that  introduces  a 
velocity  component.  This  analysis  starts  with  two  equations  listed  below:  the  continuity  equation 
for  electrons  (the  same  for  ions  assuming  charge  neutrality  and  only  a  single  ion  species:  O"^)  and 
the  phase  path  derivative  (Doppler  fi’equency). 


Where  the  phase  path  P  is: 

dt 


+  S7-in^y)  =  Q-L 


5f  =  --— 
^dt 


where: 


p(t)= 


zr(t) 

J  p(z,t)dz 


0 


80.6  n, 

f' 

•^0 


(4) 


(5) 


(6) 


and  Zr  is  the  reflection  height. 


Combining  these  two  equations  for  the  vertical  component  of  the  motion  leads  to  the  equation 
for  the  vertical  plasma  drift  velocity: 


V,  =■ 


2f 


']UQ-L)dz 


J  //  rir 


M  & 


where : 


(7) 
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in  terms  of  the  sounder  drift  velocity  v,/,  the  solar  ionization  production  rate  Q  and  the  electron 
loss  rate  L. 

This  is  an  important  equation  relating  the  vertical  plasma  drift  velocity  with  the  sounder 
measured,  Doppler  derived,  vertical  velocity  measurement.  Analysis  of  this  equation  is  plarmed 
for  future  efforts.  Using  models  of  solar  flux  as  a  function  of  wavelength  and  ionization  loss 
models  it  will  be  possible  to  compare  the  derived  plasma  drift  with  the  incoherent  scatter  radar 
measurements. 

3.  PERFORMANCE  OF  HF  COMMUNICATION  SYSTEMS 

3.1  Introduction 

Two  programs,  OpSEND  and  UPOS,  have  been  designed  to  predict  the  expected  HF 
performance  of  operational  HF  communications  and  radar  systems,  respectively.  As  inputs,  the 
two  programs  use  the  PRISM  ionospheric  model  that  is  updated  at  any  time  based  on  real-time 
measurements  fi'om  a  variety  of  sensors.  These  include  ionospheric  soxmders,  total  electron 
content  measmements  and  DMSP  satellite  measurements  of  particle  fluxes  that  specifically 
affect  the  character  of  the  auroral  ionosphere.  Based  on  our  earlier  work,  reported  extensively  in 
the  Annual  Report  2000  that  involved  the  development  of  both  programs,  OpSEND  and  UPOS 
model  the  performance  of  HF  communications  systems  at  any  location  by  generating  wide-area 
signal-to-noise  maps  that  operators  can  use  to  assess  their  ability  to  illuminate  a  remote  site  at 
any  particular  fi'equency.  These  performance  calculations  combine  a  rapid  2-D  ray-tracing 
technique  called  SMART,  developed  by  DERA  in  the  United  Kingdom  with  the  USAF  PRISM 
ionospheric  model. 

3.2  OpSEND  and  UPOS 

There  are  differences  between  the  two  programs.  OpSEND  is  a  one-way  communication 
simulation  programmed  to  hypothetical  receivers  located  anywhere  on  a  latitude/longitude  grid, 
in  all  directions  around  a  specified  transmitter  site.  These  programs  include  all  losses  including 
geometrical  spreading,  D-region  absorption  as  well  as  a  firequency/time  dependent  noise  model 
(atmospheric,  manmade  and  galactic).  The  program  computes  signal  strength,  and  signal-to- 
noise  ratio  at  each  grid  point  location.  The  second  program,  UPOS,  simulates  the  performance 
of  an  HF  backscatter  radar  system.  This  requires  a  two-way  path  fi-om  the  radar  to  the  remote 
ground  cell  and  then  scatter  back  to  the  radar.  This  requires  an  additional  parameter  for  the 
ground  backscatter  coefficient.  The  UPOS  program,  as  developed  here,  involved  only  the 
ground  backscatter  component,  although  the  actual  HF  radar  (SuperDARN)  that  was  used  to 
compare  results  was  primarily  employed  as  a  detector  of  backscatter  fi’om  field-aligned 
ionospheric  irregularities. 

The  OpSEND  program  was  delivered  to  USAF  and  installed  for  operational  use  in  the 
year  2000.  Although  the  OpSEND  product  appeared  to  work  well,  there  was  no  independent 
check  on  the  effectiveness  of  the  performance  prediction.  On  the  other  hand,  with  the  UPOS 
program,  it  was  possible  to  validate  the  approach  used  in  these  simulations  by  comparing 
performance  predictions  with  actual  radar  measurements.  The  difficulty  imposed  by  the  UPOS 
program  was  that  high  latitude  ionosphere  was  involved  in  these  comparative  tests.  The 
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SuperDARN  radar  operated  by  the  University  of  Alaska  on  Kodiak  Island  was  the  available  high 
latitude  HF  radar  system  used  here  for  the  comparative  measurements.  This  radar  scans 
continuously  over  Alaska  looking  at  both  ionospheric  scatter  from  irregularities  in  the  E  and  F- 
layers  and  ground  scatter  after  the  radar  energy  reflects  from  the  ionosphere  and  backscatters 
from  the  “rough”  ground.  For  the  UPOS  study,  it  was  the  groimd  scatter  that  the  program 
attempts  to  simulate. 

USAF  PRISM  provided  the  ionospheric  models  used  in  all  the  simulations,  both 
operational  and  for  validation.  The  PRISM  model  provides  a  latitude/longitude  grid  of  electron 
density  profiles  over  the  globe  where  the  appropriate  section  covering  the  desired  area  for 
OpSEND  is  selected  by  the  operator  and  over  the  Alaskan  region  for  UPOS.  The  results  of  these 
efforts  are  described  in  considerable  detail  in  the  Annual  Report  2000  and  only  this  brief  review 
is  included  here  although  considerable  effort  was  expended  on  these  projects.  A  typical  example 
of  the  OpSEND  operational  performance  map  is  included  here  to  illustrate  the  scope  of  the 
effort. 
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Figure  10.  Typical  HF  Illumination  Map  —  for  a  transmitter  located  near  Ottawa,  Canada 
transmitting  at  7.335  MHz  with  a  10  kW  transmitter.  The  darker  shades  of  gray  represent 

increasing  received  signal  power  (dBW). 


This  map  was  generated  for  20  Sept.  2002  at  1900  UT,  which  is  daytime  over  the  area 
shown.  The  signal  falls  off  with  distance  from  the  transmitter  indicating  a  useable  range  of  about 
1500  km  in  the  da5dime.  This  service  range  would  increase  at  night  with  a  skip  zone  surrounding 
the  transmitter.  The  OpSEND  and  UPOS  programs  handle  all  this. 
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4.  CALIBRATIONA^ALIDATION  OF  SATELLITE  BORNE  UV  ELECTRON 
DENSITY  MEASUREMENTS 

Two  far-ultraviolet  (FUV)  remote  sensing  instruments  are  planned  to  fly  on  future  DMSP 
5D3  (Defense  Meteorological  Satellite  Program)  satellites.  One  is  the  Special  Sensor  Ultraviolet 
Spectrographic  Imager,  or  SSUSI  (http://www.cpi.coni/projects/ssusi.html),  which  will  make 
disk  and  limb  observations  of  auroral  emission,  dayglow,  and  nightglow.  The  other  is  the  Special 
Sensor  Ultraviolet  Limb  Imager,  or  SSULI,  (http://www.pxi.com/praxis_publicpages/ 
SSULLhtml)  designed  to  measure  the  natural  airglow  radiation  from  atoms,  molecules,  and  ions 
in  the  upper  atmosphere  and  ionosphere.  Ionospheric  electron  density  profiles  in  the  auroral  E 
region  and  in  the  F  region  are  expected  data  products. 

Calibration  and/or  validation  of  this  space  borne  UV  technique  with  established  and 
proven  techniques  are  required.  It  was  decided  to  use  the  electron  density  profiles  available  from 
the  worldwide  network  of  digisondes  (Figure  1 1)  for  this  calibration/validation  effort. 
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Figure  11.  The  digisonde  network 


4.1  Digital  lonogram  Data  Base  (DIDBase) 

A  new  information  system  has  been  developed  that  actively  manages  data  from  the 
distributed  digital  ionosonde  network  in  support  of  a  variety  of  ionospheric  modeling  and 
comparison  projects.  As  a  result  of  the  progress  in  interactive  data  processing  technologies, 
database  management,  remote  control,  and  networking  solutions,  for  the  first  time  a  working 
enviroiunent  has  been  created  that  accesses  the  wealth  of  automatically  scaled  and  manually 
validated  ionospheric  data  derived  from  ionograms.  With  the  digisonde  data  roster  in  mind, 
DIDBase  is  expandable  to  support  other  ionosonde  models  that  use  the  universal  Standard 
Archival  Output  (SAO)  format  for  ionospheric  characteristics,  and  the  non-specific  Binary  Large 
OBject  (BLOB)  technology  to  store  the  source  ionogram  data.  Inspired  by  the  SPDDR  (Space 
Physics  Interactive  Data  Resource)  Web  project,  developed  and  maintained  by  the  WDC-A,  the 
DIDBase  expands  the  database  features  to  allow  storage  of  the  binary  ionograms  together  with 
multiple  scaling  records  for  each  ionogram.  By  relaxing  the  conventional  scenario  “one 
ionogram  -  one  set  of  derived  characteristics”,  the  DIDBase  creates  new  possibilities  for 
multiple  trace  identifications,  storage  of  alternative  ionogram  interpretations,  and  a  variety  of 
comparative  studies  of  the  automated  ionogram  processing  quality.  In  particular,  the  new 
technology  made  it  possible  to  store  several  simultaneously  observed  auroral  E  layers,  which  is 
of  special  interest  to  high  latitude  ionospheric  research.  To  provide  end  users  with  fast  access  to 
the  best  ionogram  interpretation  results,  the  DIDBase  reduces  all  versions  of  ionogram  scalings 
to  a  single  subset  of  characteristics  using  a  hierarchy  of  expert  ratings  and  quality  flags.  The  new 
information  system  includes  the  “Digital  lonogram  Data  Base”  (DIDBase),  the  expert  ionogram 
scaling  software  tool  “SAO  Explorer”,  and  the  “Automated  Data  Request  Execution  Subsystem” 
(ADRES)  to  manage  specialized  data  requests  and  monitor  their  status  from  the  inception  to  the 
final  report. 

To  compare  ionospheric  measurements  made  from  the  groimd  with  satellite 
measurements,  it  is  important  to  enhance  the  cadence  of  ground  observations  during  times  when 
the  sub  satellite  point  is  close  to  the  ground  station.  During  routine  operation,  many  of  the 
digisondes  operate  at  a  15  min  or  30  min  cadence,  some  at  60  min.  During  satellite  passes  over  a 
digisonde  station,  the  cadence  must  increase  to  one  ionogram  every  5  min  for  about  one  hour. 
Based  on  the  satellite  orbit  predictions,  the  operating  schedule  of  the  38  digisondes  around  the 
globe  participating  in  the  “CAI7VAL”  effort  must  be  adjusted  to  the  higher  cadence  for  the 
duration  of  the  passes,  in  cooperation  with  the  station  owners/operators.  Table  1  shows  the  list 
of  the  currently  activated  CALA^AL  stations.  We  have  obtained  permission  to  use  the  data  from 
these  stations  and  to  change  their  modes  of  operation. 
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Table  1.  DIDBase  Active  Digisonde  Stations 

DIDB  Active  Station  List  for  2002 


URSI 

code 

HSME 

1. 

ASOOQ 

ASCENSION  ISLAND 

2. 

RTIBB 

ATHENS 

3. 

BV530 

BUNDOORA 

4. 

CflJ2M 

CACHOEIRA  PADLXSTA 

5. 

EtL052 

CHILTON 

6. 

C0764 

COLLEGE  AK 

1 . 

DB049 

DODBBES 

8. 

DS932 

DYESS  AFB 

9. 

EG931 

EGLXN  AFB 

10. 

EA036 

EL  ABENOSHiLO 

11. 

FF051 

FAIBFORD,  UK 

12. 

FZAOM 

FORTAr^ZA 

13. 

G8J53 

GOOSE  BAY 

14. 

GR13L 

GBAHAMSTOMN 

15. 

HM19 

HAINAN 

16. 

JI91J 

JICAMABCA 

17. 

JR055 

JOLXUSKUU 

UBSX  code 

NHdE 

18. 

KS7S9 

KING  8AIMDN 

19. 

LV12P 

LCMJISVALE 

20. 

MBJ45 

MILLSTONE  HILL 

21. 

OK426 

QKINAHA 

22. 

SN437 

OSAN  AB 

23. 

PSJ5J 

PORT  STANLEY 

24. 

PA836 

PT  ABGUELLO 

25. 

TBJ77 

OAANAAQ 

26. 

PBJ18 

BCUE 

27. 

EB040 

ROQBETES 

28. 

VT139 

SAN  VTTO 

29. 

SAAOK 

SiY)  LOIS 

30. 

SMJ67 

SQNIXtESTROM 

31. 

TR169 

TE0K5O 

32. 

TO937 

NALIiOPS  IS 

33. 

WH430 

miHAN 

Also  available:  Anyang,  Gakona,  Kokubunji,  Tucuman,  Ramey 


A  total  of  approximately  60,000  ionograms  is  expected  for  a  one-month  CAIW AL  period.  This 
data  volume  requires  the  development  of  a  “digital  ionogram  database  or  DIDBase,  with  direct 
access,”  to  SAO-Explorer,  the  fourth  generation  tool  for  digisonde  data  visualization  SAO-X. 
DIDBase  was  developed  in  consultation  with  Terence  Bullett  of  AFRL,  and  is  briefly  described 
below.  A  second  requirement  is  a  mechanism  that  automatically  reads  incoming  requests  from 
the  satellite  operators  and  adjusts  programs  and  schedules  of  the  Internet-enabled  digisondes. 
SAO-X  was  upgraded  to  include  die  ADRES  subsystem  to  handle  this  task. 
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Figure  12.  Organization  of  DIDBase 


Figure  12  shows  the  organization  of  the  data  flow  for  the  CALA^AL  operation.  The  digisonde 
data  are  deposited  in  the  database,  and  the  SAO-X  workstations  have  direct  access  to  these  data. 
This  is  important  since  the  ARTIST  autoscaling  of  ionograms  [Reinisch  and  Huang,  1983]  is  not 
always  error  free,  and  manual  editing  of  the  scaling  is  required  to  assure  correct  electron  density 
profiles  for  comparison  with  the  UV  results.  Two  examples  of  autoscaled  ionograms  from 
Jicamarca  for  5  Dec  2002  are  shown  in  the  left  column  in  Figure  13.  The  autoscaling  at  2100  UT 
(upper  panel)  is  very  accurate  with  values  of  14.2  MHz  and  419  km  for  foF2  and  hmF2.  Manual 
editing  (right  column)  changed  these  values  by  only  ~1%,  although  the  calculated  vertical  total 
electron  content  (TEC)  changed  by  5%.  At  2315  UT  ARTIST  was  misled  by  a  small  echo  gap  in 
the  trace  (possibly  caused  by  interference)  and  grossly  underestimated  foF2,  hmF2,  and  TEC,  as 
can  be  seen  in  the  edited  ionogram  on  the  right.  While  ARTIST  has  been  slightly  improved  since 
this  record  was  made  to  avoid  hanging  up  on  small  gaps,  these  mistakes  are  still  occurrmg  for  the 
DISS  stations  where  the  AF  Frequency  License  creates  wide  gaps  in  the  traces. 


Figure  13.  Illustrating  the  need  for  manual  editing 


The  function  of  the  ADRES  subsystem  is  illustrated  in  Figure  14.  ADRES  consists  of  two  Java 
based  programs,  DIDBFill  (loading  data  into  DIDBase)  and  DIDBReqPro  (DIDBase  Request 
Processor).  The  ADRES  subsystem  can  request  validated  ionospheric  data  in  form  of  plain  text 
reports.  If  the  requested  data  are  already  in  DIDBase  and  had  been  validated,  the  ADRES 
generates  the  report  immediately.  A  provision  is  made  to  manage  the  requests  for  special  modes 
of  ionosonde  operation  (e.g.,  high  ionogram  rate  during  a  satellite  pass  over  the  station,  a 
coordinated  campaign,  or  an  event  of  interest)  and  data  that  are  not  available  in  DIDBase  or  are 
not  manually  validated.  The  ADRES  has  a  mechanism  to  automatically  read  incoming  requests 
to  adjust  programs  and  schedules  of  the  Internet-enabled  Digisondes.  To  acquire  Digisonde  data, 
the  ADRES  maintains  a  list  of  FTP  servers  where  the  ionogram  data  can  be  found.  A  number  of 
digisonde  stations  deliver  their  real  time  data  directly  to  a  WDC  and  to  DIDBase.  As  soon  as  the 
data  are  ingested,  a  message  is  generated  to  the  SAO-X  operators  to  validate/edit  the  autoscaling 
results.  When  the  quality  control  procedure  is  completed,  the  final  report  is  generated  and 
delivered  to  the  requesting  party. 
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Figure  14.  The  functions  of  the  DIDBase  ADRES  subsystem 
4.2  CALA^AL  Dry  Run  with  TIMED/GUVI 


The  TIMED  spacecraft  carries  a  UV  instrument,  GUVI,  which  is  almost  identical  to  the 
SSUSI  instrument  to  be  flown  on  DMSP.  A  real  time  dry  run  was  initiated  for  the  TIMED 
spacecraft  orbit  during  a  3 -day  campaign  fi'om  25-27  April  2002  to  test  and  exercise  the 
DIDBase/digisonde  network  operation.  Ionospheric  model  predictions  for  this  magnetic  storm 
period  were  compared  with  the  characteristics  measured  by  the  digisondes  at  1 1  stations  in  the 
Western  hemisphere  (Table  2). 

Table  2.  Eleven  digisonde  stations  were  included  in  the  TIMED  dry  run 


# 

URSI  CODE 

STATION  NAME 

1 

WP937 

Wallops  Island 

2 

DS932 

Dyess 

3 

ASOOQ 

Ascension  Island 

4 

MHJ45 

Millstone  Hill 

5 

JI91J 

Jicamarca 

6 

PSJ5J 

Port  Stanley 

7 

SAAOK 

SAO  LUIS 

8 

CAJ2M 

CACHOEIRA  PAULI STA 

9 

FZAOM 

Fortaleza 

10 

EG931 

Eglin 

11 

PAS  3  6 

Pt  Arguello 

Model  simulations  (gray)  of  the  diurnal  variations  of  the  F2  layer  peak  density  and  height  are 
shown  in  Figure  15  for  six  digisonde  locations.  The  measurements  are  shown  as  black  curves. 
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Figure  15.  Simulations  (gray)  of  NmF2  (cm'^  x  10'®)  on  the  left  and  hniF2  (km),  on  the  right,  for 
the  Southern  hemisphere,  are  compared  with  digisonde  measurements  (black). 

(courtesy  of  Paul  Straus,  Aerospace  Corporation) 

The  digisonde  software  (both  DPS-4  and  DGS-256)  has  been  modified  so  that  it 
automatically  responds  to  CALA^AL  campaign  requests.  Requests  received  fi-om  Aerospace  are 
entered  into  the  ADRES  system  to  automatically  reprogram  the  sounding  schedule  and 
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measurement  programs  of  selected  stations.  A  request  example  for  25  April  2002  is  shown  in 
Table  3  (all  times  are  given  in  UT). 


Tables.  Request  for  CALA^AL  campaign 


ADD 

WP937/ 

DS932, 

ASOOQ, 

MHJ45, 

JI91J, 

PSJ5J, 

PSJ5J, 

SAAOK, 

CAJ2M, 

FZAOM, 


2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 


21:09:00, 

00:05:00, 

17:43:00, 

21:11:00, 

20:56:00, 

17:28:00, 

19:08:00, 

19:21:00, 

19:16:00, 

19:21:00, 


2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 

2002.04.25 


22:09:00 

01:05:00 

18:43:00 

22:11:00 

21:56:00 

18:28:00 

20:08:00 

20:21:00 

20:16:00 

20:21:00 


ADRES  also  issues  text  reports  to  Aerospace  containing  the  most  relevant  measurements 
(an  example  is  given  in  Table  4).  Starting  in  late  2002,  these  reports  used  the  standard  SAO-4 
data  format  (http://129.63. 134.1/fi*amesd.htm) 


Table  4.  CAL/VAL  campaign  request  for  the  Millstone  Hill  digisonde 

UMass  Lowell  Ionospheric  REPORT  REQUESTED  BY  DMSP  Space  Science 

Applications  Laboratory  The  Aerospace  Corporation 

FOR  STATION:  MILLSTONE  HILL,  URSI:  MHJ45 

START  TIME:  2002“Apr-26  (116)  21:29:00  UT 

END  TIME:  2002-Apr-26  (116)  22:29:00  UT 


Auto  scaled 

-  A, 

Validated  -  V, 

TEC  units  - 

[10^16  m 

^-2] 

yyyy . MM . dd 

(DDD) 

HH:inm:ss 

A/V 

foF2 

[MHz] 

]iinF2 

[km] 

TEC 

[TECU] 

foEa 

[MHz] 

2002.04.26 

(116) 

21:30:00 

V 

9.750 

308.0 

32.5 

— 

2002.04.26 

(116) 

21:35:00 

V 

9.750 

307.6 

32.2 

— 

2002.04.26 

(116) 

21:40:00 

V 

9.750 

307.6 

31.6 

— 

2002.04.26 

(116) 

21:45:00 

V 

9.750 

307.8 

31.2 

— 

2002.04.26 

(116) 

21:50:01 

V 

9.850 

314.3 

33.3 

— 

2002.04.26 

(116) 

21:55:00 

V 

9.850 

311.5 

32.2 

— 

2002.04.26 

(116) 

22:00:00 

V 

9.850 

309.9 

31.2 

— 

2002.04.26 

(116) 

22:05:00 

V 

9.950 

312.7 

32.6 

— 

2002.04.26 

(116) 

22:10:00 

V 

9.950 

313.2 

32.7 

— 

2002.04.26 

(116) 

22:15:00 

V 

9.850 

307.9 

30.0 

— 

2002.04.26 

(116) 

22:20:00 

V 

9.950 

312.7 

31.9 

— 

2002.04.26 

(116) 

22:25:00 

V 

9.950 

309.0 

30.0 

— 
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4.3  Auroral  E  Layer  Studies  at  Sondrestrom 

It  is  planned  to  deduce  the  E  region  electron  density  profiles  from  the  ni^ttinie  SSUSI 
measurements.  Of  special  interest  will  be  the  ionization  produced  by  electron  precipitation  in  the 
auroral  zone.  This  ionization  also  shows  on  ionograms  in  the  form  of  auroral  E  layers,  as 
illustrated  in  Figure  16.  Enhanced  retardations  of  the  echoes  occur  at  1.7  MHz,  5.2  MHz,  and  5.8 
MHz,  indicating  ionization  peaks  within  the  field  of  view  of  the  digisonde’s  antenna  system. 
With  a  beam  width  of  approximately  1  radian,  an  area  of  50  km  radius  at  E  region  heights 
contributes  to  the  ionograms.  This  area  contains  regions  of  higher  density,  3.3x10  and  4.2x10 
m■^  and  the  background  ionization  of  3.6xl0'°  m'l  At  times  when  only  a  single  auroral-E  layer 
is  present,  it  is  likely  to  be  associated  with  diffuse  aurora.  The  footprint  of  the  SSUSI  instrument 
is  about  a  factor  of  2  smaller  than  the  digisonde’s,  but  it  still  will  contain  significantly  varying 
d.6iisitics 

SAO  Explorer,  V3.2.03b5  M00.09.26  (270>  21:05:00  SL 
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Figure  16.  E  region  ionograms  at  Sondrestrom  on  26  Sep  2000  at  2105  UT.  The  data  pps  at  2.9 
MHz  and  5.4  MHz  are  caused  by  transmission  restrictions  in  the  AF  Frequency  License 

The  collocated  incoherent  scatter  radar  at  Sondrestrom  was  making  profile  measurements 
during  this  time  in  a  fixed  beam  direction  parallel  B.  The  ISR  beam  width  is  0.6°  corresponding 
to  a  1  km  footprint  in  the  E  region.  The  one-minute  sequence  of  ISR  profiles  displayed  in  Figure 
17  shows  an  order  of  magnitude  variation  in  NmE.  The  digisonde  measurements  of  NmE  firom  the 
ionogram  at  2105  UT  were  added  to  the  figure.  While  the  digisonde  “sees”  all  densities  at  once, 
the  ISR  sees  the  blobs  of  ionization  move  through  its  beam  as  function  of  time. 
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Figure  17.  One-minute  sequence  of  incoherent  scatter  radar  electron  density  profiles  (black)  at 
Sondrestrom.  One  digisonde  profile  fi-om  the  2105  UT  ionograms  is  shown  in  red,  superimposed 
on  the  2105  UT  ISR  profile.  Two  other  peak  densities  (red  marker)  from  the  2105  UT  ionograms 
are  shown  on  the  profiles  with  similar  NmE.  (ISR  profiles  courtesy  Jeff  Thayer,  SRI). 
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6.  SUPPORT 


During  the  lifetime  of  this  contract,  the  Center  for  Atmospheric  Research  (CAR)  has 
provided  extensive  support  to  the  USAF  and  its  facilities  aroimd  the  world.  Below  is  a  list  of  the 
travel  carried  out  by  CAR  as  part  of  this  support  function. 


9-13  Feb  1996 

I- 6  June  1996 

II- 17  Nov  1996 
17-21  Mar.  1997 


22-29  Sept.  1997 
21-22  Jan  1998 
19  Feb-4  Mar  1998 
21  Mar  -  4  Apr.  1998 


28  Apr -2  May  1998 

14- 21  Aug  1998 

28  Aug  1998 
18  Sept.  1998 
23-25  Sept.  1998 
28-29  Oct.  1998 

15- 20  Nov.  1998 

29  Nov -2  Dec.  1998 
29  Nov -2  Dec  1998 
28-31  Dec  1998 

4-8  Jan  1999 
27-31  Dec  1999 
27-31  March  2000 
18-24  June  2000 
9-12  Oct.  2000 


Bermuda 
Osan  AB,  Korea 
Tromso,  Norway 
RAF  Fairford,  UK 
USAF  AWS/Space  Forecast 
CLRC  RAL 

Sondrestrom,  Greenland 

Wallops  Island,  VA 

McClellan,  Vandenburg  AFB,  CA 

Korea 

Japan 

Eglin  AFB,  FL 

Dyess  AFB,  TX 

Goose  Bay,  Canada 

Narsarsuaq,  Greenland 

Wallops  Island,  VA 

Goose  Bay,  Canada 

Wallops  Island,  VA 

Wallops  Island,  VA 

Sacramento,  CA 

Eglin  AFB-Ft.  Walton  Beach,  FL 

Ramey  Solar  Observatory,  San  Juan,  PR 

Bermuda 

Sheep  Creek  and  Gakona,  AK 

Bermuda 

Tobyhanna,  PA 

Tobyhanna,  PA 

Bermuda 
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11-18  Mar.  2001 
4-6  March  2002 
7-12  May  2002 
3-10  July  2002 
15-22  Sept.  2002 


Ascension  Island 
Wallops  Island,  VA 
HAARP  -  Gakona,  AK 
Tucuman,  Argentina 
Campo  Grande,  Brazil 
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